Deconfining Phase Transition on Lattices with Boundaries at Low Temperature 
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In lattice gauge theory (LGT) equilibrium simulations of QCD are usually performed with periodic 
boundary conditions (BCs). In contrast to that deconfined regions created in heavy ion collisions 
are bordered by the confined phase. Here we discuss BCs in LGT, which model a cold exterior of 
the lattice volume. Subsequently we perform Monte Carlo (MC) simulations of pure SU(3) LGT 
with a thus inspired simple change of BCs using volumes of a size comparable to those typically 
encountered in the BNL relativistic heavy ion collider (RHIC) experiment. Corrections to the usual 
LGT results survive in the finite volume continuum limit and we estimate them as function of the 
volume size. In magnitude they are found comparable to those of including quarks. As observables 
we use a pseudocritical temperature, which rises opposite to the effect of quarks, and the width of 
the transition, which broadens similar to the effect of quarks. 

PACS numbers: PACS: 05.10.Ln, 11.15.Ha 



I. INTRODUCTION 

At a sufficiently high temperature QCD is known to 
undergo a phase transition from our everyday phase, 
where quarks and gluons are confined, to a deconfined 
quark-gluon plasma. Since the early days of lattice gauge 
theory simulations of this transition have been a subject 
of the field flj, see for reviews. Naturally, such sim- 
ulations focused on boundary conditions (BCs), which 
are favorable for reaching the infinite volume quantum 
3 continuum limit quickly. On lattices of size A',- A^^ 
these are periodic BCs in the spatial volume V = {a Ng)^, 
where a is the lattice spacing. For a textbook, see, e.g., 
Ref. [3. 

The physical temperature of the system on a A^r A^^ 
lattice, Nr < Ns, is given by 



T - 



aNr 



1 



(1) 



where a is the lattice spacing. In this paper we set the 
physical scale by 



= 174 MeV 



(2) 



for the deconfinement temperature, which is approxi- 
mately the average from QCD estimates with two light 
flavor quarks [2] in the infinite volume extrapolation. The 
relation ([T]) implies for the temporal extension of the sys- 
tem 



Lr = aN-r = 1.13 fermi . 



(3) 



For the deconfinement phase created in a RHIC the in- 
finite volume limit Ns/Nj- —f 00 for fixed Nr and sub- 
sequently Nt — > 00 (Lr — a Nr finite) does not apply. 
Instead we have to take the continuum limit as 



and periodic BCs are incorrect because the outside is 
in the confined phase at low temperature. Details are 
discussed in the next section. 

In collisions at the BNL RHIC one expects to cre- 
ate an ensemble of differently shaped and sized volumes, 
which contain the deconfined quark-gluon plasma. The 
largest volumes are those encountered in central colli- 
sions. A rough estimate of their size is 

TT X (0.6 X Au radius)^ x c x (expansion time) 

= (55 fermi'^) x (a few fermi) (5) 

where c is the speed of light. To imitate this geometry, 
one may want to model spatial volumes of cylindrical and 
other geometries. In our exploratory study we do not try 
to imitate a realistic ensemble of deconfined volumes, but 
are content to estimate the magnitude of corrections one 
may expect. So we stay with N^ volumes and focus on 
results in the continuum limit for 



L, = aNs = (5 - 10) fermi 



(6) 



Ng/Nr = finite , A^,- ^ 00 , Lr finite , 



(4) 



Finite volume corrections to the infinite volume contin- 
uum limit are expected to be relevant as long as the vol- 
ume is not large compared to a typical hadronic correla- 
tion length, which is about one fermi. For relatively small 
volumes an appropriate modeling of the BCs is necessary. 

In the next section we introduce BCs, which reflect a 
(very) low temperature outside of the deconfined region. 
Two constructions, the "disorder wall" and the "confine- 
ment wall" are carried out and shown to exhibit the usual 
asymptotic scaling properties in the finite volume con- 
tinuum limit. While the confined wall is physically more 
accurate, the disorder wall can be easily implemented in 
MC calculations. For the latter MC calculations of pure 
SU(3) LGT are performed in sect ion [Hi] and evidence for 
scaling is already found for systems sizes which can be 



2 



simulated on PC clusters. Summary and conclusion are 
given in section [TVl 



II. BOUNDARY CONDITIONS 

Statistical properties of a quantum system with Hamil- 
tonian H in a continuum volume V, which is in equi- 
librium with a heatbath at physical temperature T, are 
determined by the partition function [4| 



Z{T,V)=Tve-"/^ = J2{ 



-H/T\ 



(7) 



where the sum extends over all possible states \4>) of the 
system and the Boltzmann constant is set to one. Impos- 
ing periodic boundary conditions in Euclidean time r and 
bounds of integration from to l/T, one can rewrite the 
partition function ([7]) in the path integral representation: 



Z{T,V) = Dcjjcxp^- 



l/T 



cLtL} 



(8) 



Nothing in this formulation requires to carry out the in- 
finite volume limit. In the contrary, if one deals with 
rather small volumes for which fluctuations of the mean 
values are not negligible, the idea to consider V ^ oo 
instead of V appears to be rather obscure. 

An obvious problem of applying equilibrium thermo- 
dynamics to deconfined volumes at RHIC is that there 
is no heatbath in sight with which the system could be 
in equilibrium. However, arguments have been made in 
the literature that after the rapid heating quench, when 
the deconfined volume is at about its maximum size, a 
(pseudo) equilibrium state is reached for a transitional 
period, which is reasonably long on the scale of the re- 
laxation times involved. These arguments are not be- 
yond doubts (some are raised in one of our own papers 
[7|), but they are not a topic of our present work. Here 
our assumption is that there is some truth to the belief 
that the bulk properties can be described by equilibrated 
finite temperature QCD. 

In MC simulation the updating process provides the 
heatbath and finite volumes can be equilibrated with all 
kind of BCs imposed, although in practice most simula- 
tions have used periodic BCs. For simplicity and to be 
definite we restrict our discussion to pure SU(3) LGT. 
Generalization of the arguments of this section to full 
QCD appears to be straightforward. We use the Wilson 
action given by 



S{{U}) 



(3s 



^ReTr(t/n) 



(9) 



[/□ — Ui-^j-^Ujj^i^Ui^j^Uj^i-^, where the sum is over all pla- 
quettes of a 4D simple hypercubic lattice, ii, ji, 12 and 
j2 label the sites circulating about the plaquette and Uji 
is the S'[/(3) matrix associated with the link {ij). The 



reversed link is associated with the inverse matrix, and 
13^ is related to the bare coupling constant by (3^ = 6/5^. 
The theory is defined by the expectation values of its 
operators with respect to the Euclidean path integral 



(10) 



where the integrations are over the invariant group mea- 
sure, which was for compact groups like SU(3) introduced 
by Hurwitz 9] (Haar generalized it later). 

Numerical evidence suggests that for Nt fixed and 
A^s — > 00 SU(3) lattice gauge theory exhibits a decon- 
fining phase transition (sub- or superscript t stand for 
transition) at some coupling (31{Nt) — 6/g'^{Nr), which 
is weakly first order [Tl|. The scaling behavior of the 
deconfining temperature is 



(11) 



where ct is a constant and we use the lambda lattice 
scale 

aAi = = Kg') (605')"'^^^"°' e-V(2fco/) ^ 

(12) 

where a is the lattice spacing. The coefficients bo and 
61 are perturbatively determined by the renormalization 
group equation and independent of the renormalization 
prescription [l3 |. 



11 3 , , 34/3 

oq = — ^ „ „ and oi = 



3 167r2 



3 Vl67r2 



(13) 



For perturbative and non-perturbative corrections we 
adopt the analysis of fT3] in the parameterization of Q: 



1 + aie 



-0.2/9 



0.3 g +043 



(14) 



with ai = 71553750, 02 = 19.48099, 03 = -0.03772473, 
and 04 = 0.5089052. 

We want to model an equilibrium situation surrounded 
by cold boundaries at 300 K. The effects are expected 
to penetrate at least a few correlation lengths into the 
volume. Note that this is the correlation length set by 
hadronic interaction, which can be defined by an appro- 
priate inverse mass. On the lattice ^/a governs the con- 
tinuum limit. This should not be confused with corre- 
lations due to the phase transition. As the lattice reg- 
ularization allows to construct finite volume continuum 
limits as well as the infinite volume continuum limit, the 
question is: Which construction has the best chances to 
capture the physics realistically for a quark-gluon plasma 
in a small volume as typically created at RHIC? 

Let us first exemplify that LGT thermodynamics al- 
lows not only to approach the infinite volume continuum 
limit, but also finite volume continuum limits. For in- 
stance, we can define the thermodynamics on a torus 
which has the volume of, say, (lOfermi)'^. To achieve 
a continuum limit, we have to send the lattice spacing 



3 



a ^ in units of the physical scale. This is governed 
by the renormalization group equation and requires in- 
finitely many lattice points Ns — Ls/a —>■ oo, while the 
physical volume of the lattice stays finite by arranging 
Ls/fermi — ci, where ci is a constant, e.g., ci = 10. The 
temperature of such a system is regulated by choosing 
Lr/Ls = C2, where C2 is a second constant. For describ- 
ing deconfined volumes at RHIC the thus defined toroidal 
mini-universe is even less suitable than the conventionally 
used infinite volume limit, because correlations are artifi- 
cially propagated through the periodic BCs. We present 
now two constructions of BCs, which reflect cold exterior 
volumes. The second is physically more realistic, but nu- 
merically more difficult to implement. 



The disorder wall allows for a finite volume continuum 
limit, which is approached for /3g — > oo, and the usual 
scaling law holds. Eventually (obviously for /3q > 24) the 
value of (3f increases also, unless the outside temperature 
is exactly zero and, hence, f3f — 0. An outside temper- 
ature of 300 K is on the scale of the inside temperature 
so close to zero that in practical MC simulations the ef- 
fects of the (3f increase are not noticeable and /3f = is 
a safe choice independently of Pq. Simulation results in 
section HII] support that scaling in /3q holds already for 
/3q values, which can be reached in practice. 



B. Confinement wall 



A. Disorder wall 

Imagine an almost infinite space volume V = L^, which 
may have periodic BCs, and a smaller (very large, but 
small compared to V) sub- volume Vq = L^q. The com- 
plement to Vb in will be called Vi . The number of tem- 
poral lattice links Nt is the same for both volumes. We 
want to find f3^ values and lattice dimensions so that scal- 
ing holds, while Vq is at a temperature of Tq = 174 MeV 
and Vi at Ti = 300 K. We set the couphng to 
for plaquettes in Vq and to /3f for plaquettes in Vi. For 
that purpose any plaquette touching a site in Vi is con- 
sidered to be in Vi. Let us take /3f = 5.7, which is 
at the beginning of the SU(3) scahng region. We have 
174 MeV = 2.02 x lO^^K and, therefore. 



To 



2.02 X 10^2 
300 



ao 



(15) 



where ai is the lattice spacing in Vi, 
the scaling equation 



0,1. Using m 



300/A(/3f) = 2.02xlOiVA(/3o') 



(16) 



yields (3^ = 



24.496. Tr estimates from MC calculations 



10 , 



of the literature extrapolate then to Lr > 2.74 x 10 
for the temporal lattice extension needed for a decon- 
fined phase at /3g = 24.5. This illustrates the orders of 
magnitude involved. 

So, in practice we can only have /Jq, but not /Jf in the 
scaling region, say /3q = 6. There is still a reasonable 
choice for (3\. The scaling argument (jTS]) shows that in 
lattice units a correlation length ^, say one fermi, is at 
/3q much larger than at /Jf: 



(e/«i) 



^«10i°. 



(17) 



Now ^/ao is of order one at /3q ~ 6, so that ^/ai becomes 
very small and the lowest order of the strong coupling ex- 
pansion applies (see 0, 13] for the string tension and 
for the glueball mass) in which ^/ai ^ — 1/ ln(/3f /18) so 
that /3f K, 10"^" . We call this construction the disorder 
wall. 



In the continuum limit the disorder wall separates re- 
gions in which a correlation length ^ takes on entirely 
different magnitudes when measured in unit of the lat- 
tice spacing (|17p . Although we are only interested in the 
physics inside the wall, a construction for which the dis- 
tance of one fermi stays in lattice units continuous across 
the boundary is clearly more physical. Along similar lines 
as before this can be achieved by using an anisotropic lat- 
tice outside of the deconfined region. Let us denote in Vi 
the spacelike links by Os, the timelike links by a^, and 
use there the Wilson action 



S{{U}) 



PI 
3 



^ReTr(C/nJ 



PI 
3 



ReTr([/nJ 



(18) 

where /3f and are the couplings of the spacelike and 
timelike plaquettes, respectively. The lambda scale of 
this action has been inv estig ated by Karsch [l^ and in 
the continuum one finds hj 



PI 

PI 



As we aim at 



10" 



-10 



(19) 



(20) 



the resulting orders of magnitude are even more astro- 
nomical than before. The sublattice Vi is again driven 
out of the scaling region, which would be reached for suf- 
ficiently large values of /^Karsch ~ vPsPt- For all practi- 
cal purposes we are driven into the strong coupling region 
and can set P^ = 0, so that the simulation of the confined 
world becomes effectively 3D. By measuring an appropri- 
ate correlation length, for instance via the string tension 
or glueball mass, we can non-perturbatively tune /3f , so 
that as — Oo holds. A MC simulation has then to include 
an outside world with a non-zero /3f . This is physically 
quite interesting, but computationally more demanding 
than using the disorder wall, for which we present MC 
simulations in the following. 
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Periodic BCs 


Disorder wall BCs 




^ meas 




^meas 




12 


32 X 10 000 


5.6904 (27) 


64 X 20 000 


6.110(34) 


16 


32 X 10 000 


5.6912(11) 


64 X 20 000 


5.8460 (83) 


20 


32 X 10 000 


5.69184 (69) 


32 X 10 000 


5.7744 (59) 


24 


32 X 10 000 


5.69170 (41) 


32 X 10 000 


5.7426 (30) 


32 


32 X 10 000 


5.69225 (16) 


32 X 10 000 


5.7192(11) 



TABLE I: Number of measurements and pseudo-transition 
coupling estimates for 4 x N'^ lattices. 



III. SU(3) MC SIMULATIONS WITH THE 
DISORDER WALL 



In the Monte Carlo calculations of this section we ap- 
proximate a cold exterior by using the disorder wall BCs. 
In practice this means, we simply omit plaquettes, which 
involve links through the boundary. So we drop the sub- 
script in the /3q definition of the previous section and 
return to simply using the [3^ notation. For both peri- 
odic and disorder wall BCs we present an analysis of data 
from simulations on Nr x lattices with Nr — 4 and 6. 
The Ns values and our statistics in sweeps are compiled 
m tables U and m 

As in |l5 we use the maxima of the Polyakov loop 
susceptibility 



Xmax — 



1 



1^1') -(1^1)'], 



(21) 



to define pseudo-transition couplings (3p^{Ns; Nr). For 
periodic BCs, indicated by the superscript p of a^, they 
have a finite size behavior of the form 



N. 



(22) 



Fits to this form yield f3f{Nr) and estimates are given in 
Boyd et al. Our estimates for f3pf{Ns;4:) are summa- 
rized in table m Within statistical errors the lattice size 
dependence of (3pf{Ns; 4) is almost negligible. A fit of our 
data with periodic BCs to (HH) gives /3f (4) = 5.69236 (21) 
in agreement with the value 5.6925 (2) reported in p^ . 
From now on we use /3f (4) = 5.69236 (21) as the infinite 
volume limit for periodic and disorder wall BCs. 

Our estimates of /3p((iVs; 6) are summarized in table [TTl 
We took only few Nt = 6 data for periodic BCs, because 
they consume already considerable CPU time and results 
exist already in the literature. Again they show almost 
no lattice size dependence. Fitting them to ((22| gives 
/3f (6) = 5.8926 (18) in agreement with (3f{6) = 5.8941 (5) 
from [13]. In the following we use the latter, more accu- 
rate, /9f (6) estimate of the literature. 

The disorder wall BCs introduce an order iV^ distur- 





Periodic BCs 


Disorder wall BCs 




Nmeas 


Plt{Ns;6) 


N meas 


f3'r.t{Ns;6) 


18 


32 X 10 000 


5.8932 (48) 


32 X 10 000 


6.47 (14) 


20 


- 


- 


32 X 10 000 


6.27 (04) 


24 


32 X 10 000 


5.8934 (26) 


32 X 10 000 


6.089 (23) 


28 






32 X 10 000 


6.012(11) 


32 


32 X 10 000 


5.8927(12) 


32 X 10 000 


5.9812 (73) 


40 






32 X 10 000 


5.9463 (53) 


48 






16 X 12 000 


5.9271 (38) 



TABLE IL Number of measurements and pseudo-transition 
coupling estimates for 6 x Ng lattices. 
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FIG. 1: Fits of pseudo-transition coupling constant values and 
their infinite volume extrapolations. 



bance, so that Eq. 



becomes 



alti^s-.Nr) 



N, 



(23) 



N, 



N^ 

Ns 



where the superscripts d of the coefficients af indicate 
disorder wall BCs. 

In Fig. [1] we show the fit ((23|) to /3f (4) and our pseudo- 
transition values /?pj [Ns] 4) from simulations with the dis- 
order wall BCs. The very precise infinite volume estimate 
/3f (4) from simulations with periodic BCs is included in 
the disorder wall data to stabilize the fit at large volumes. 
For comparison the fit for our Nr = 4 data from sim- 
ulations with periodic BCs is also given. While finite size 
corrections are practically negligible for the simulations 
with periodic BCs, this is not the case for the disorder 
wall BCs. Q is the goodness of the fit (e.g., chapter 2.8 
of 18]). 

We also include the fit ^ to /3f (6) and our l3^t{Ns; 6) 
disorder wall data in Fig. [TJ As expected both disorder 
wall curves show strong finite lattice size effects. This 
is not automatically of physical relevance. Important is 
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FIG. 2: Estimate of finite volume corrections to the decon- 
finement temperature, set at 174 MeV for an infinite volume. 



whether universal corrections survive in the finite volume 
continuum limit. Here universal means that the correc- 
tions do not depend on the lattices used in the simula- 
tions, once these lattices are sufficiently large. In the 
following we test this for Nr = 4 and = 6 using the 
lambda scale (|12|) to calculate estimates in physical units. 
Although our (3^ values used are rather small, they are 
in the previously reported [l^l scaling region for (fO|) , so 
that it is reasonable to expect universal behavior with 
moderate corrections. 



The infinite volume T'^ value the lambda scale 
and Eq. (fTTj) give us the g"^ dependence of the lattice spac- 
ing a in units of fermi. Using the = 4 fit to of 
Fig. [T] together with Eqs. ([T]) and ^ allows to eliminate 
and we plot the resulting function T'^{Ls) in Fig. [21 
Repeating this procedure for our TVt- = 6 fit to (|23|l gives, 
as is seen in Fig. [21 almost the same T^{Ls) dependence 
and thus provides some evidence that our Nr = 4 and 6 
results are already representative for the finite volume 
continuum limit. For a box of volume (lOfermi)'^ the 
pseudocritical temperature T'^ is about 5% higher than 
the infinite volume estimate and this correction increases 
to about 17% for a (Sfermi)^ box. For comparison we 
use the same procedure for analyzing the finite size de- 
pendence obtained by fitting the Nt = 4 and 6 data with 
periodic BCs (including and enforcing the /3f (6) = 5.8941 
limit for Nr = 6) to Eq. ([221) • This gives the two lower 
curves of the figure, which fall almost on top of one an- 
other. In that case error bars are considerably smaller 
than for the disorder wall data. 

For a first order phase transition the maxima of the 
Polyakov loop susceptibility have to scale with the system 
volume to reproduce the delta function like singularity of 
a first order transition. For Nr = 4 our Xmax values are 
listed in table IIIII and for Nr = 6 in table IIVI Using 
Ns > 16 for Nr ^ i and Ns > 20 for A^^ = 6 acceptable 





Periodic BCs 


Disorder wall BCs 




Xmax 




Xmax 




12 


3.585 (71) 


0.0207 (11) 


1.715 (27) 


0.448 (18) 


16 


7.62 (16) 


0.0103 (13) 


1.879 (26) 


0.0997 (21) 


20 


16.17 (67) 


0.00498 (68) 


2.525 (84) 


0.0440 (27) 


24 


28.6 (1.1) 


0.00277 (32) 


3.58 (14) 


0.0225 (16) 


32 


73.0 (2.0) 


0.00132 (16) 


7.67 (39) 


0.00709 (61) 



TABLE III: Maxima of Polyakov loop susceptibility and width 
of the transition for lattices 4 x N?. 





Periodic BCs 


Disorder wall BCs 




Xniax 






A/32^3 


18 


2.47 (10) 


0.0303 (21) 


2.15 (20) 


0.84 (04) 


20 






1.89 (10) 


0.322 (25) 


24 


5.00 (27) 


0.0162 (15) 


1.90 (11) 


0.140 (13) 


28 






2.18 (10) 


0.0712 (69) 


32 


11.34 (55) 


0.00803 (66) 


2.82 (16) 


0.0422 (37) 


40 






4.21 (31) 


0.0234 (17) 


48 






5.44 (99) 


0.0123 (38) 



TABLE IV: Maxima of Polyakov loop susceptibility and width 
of the transition for lattices 6 x Nf. 



fits to the straight-line form 

Xmax = di+d2 [NjNr f 



(24) 



are obtained and shown together with their Q values in 
Fig. [3l (for Nr — 6 and periodic BCs only two data points 
are fitted, so that there is no Q value in that case). To 
enhance the scale for the disorder wall fits, they are dis- 
played on the right ordinate. The leading coefficients ob- 
tained from fits for disorder wall data differ from those 
for the data with periodic BCs: ^2 = 0.01190 (60) versus 
= 0.1436 (37) for A^ = 4 and = 0.00843 (95) versus 
= 0.0723 (70) for A'^,- = 6. This is possible because the 
Polyakov loop maxima are not physical observables, but 
bare quantities. 

Let us now consider the width of the transition. For 
a 4 X 16"^ lattice with disorder wall BCs the Polyakov 
loop susceptibility as a function of f3^ is shown in Fig. [H 
We use reweighting [H, [l^ to cover a range of and 
define ^Pys ^ width of a peak at 2/3 of its height. 
On large lattices the figures look less nice than Fig. [H 
but provide still sufficient information to extract Afiy^- 
Our estimates for A^,- 4 are listed in table IIIII and for 
At- = 6 in table IIVI The data are fitted to the form 



for periodic BCs and to 



Nr_ 

N, 



Nr 

N. 



Nr 

Ns 



(25) 



(26) 



6 





80 




70 




60 


O 




m 


50 


o 






g 


40 
tu 


to 




a. 


30 


X 

to 




E 


20 




10 








N^=4 periodic BCs 
N^=6 periodic BCs 
N^=4 disorder wall 
N^=6 disorder wall 




100 200 300 400 500 

(Ns/N,)2 



g 

03 
T3 

o 



E 



CO. 



0.1 
0.08 
0.06 
0.04 
0.02 



Disorder wall BCs 
Periodic BCs 



Q=0.67 



Q=0.85 





0.12 0.14 0.16 0.18 0.2 0.22 0.24 0.26 

N,/N3 



FIG. 3: Maxima of Polyakov loop susceptibility. 



FIG. 5: Fits of the iV^ = 4 width of the transition. 
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FIG. 4: The Polyakov loop susceptibility on a 4 x 16^ lattice 
with disorder wall BCs. 



for the disorder wall BCs. The first term reflects in both 
cases the delta function singularity of a first order phase 
transition. The leading order corrections to that differ 
due to the influence of the BCs. 

For Nr = 4 the final fits (see below) are shown in 
Fig. [5l From the disorder wall data we have omitted our 
smallest 4 x 12"^ lattice from the fit, because the width be- 
comes for it so broad that it spoils Q (larger lattices than 
with periodic BCs are needed). The leading order coeffi- 
cients are then cf = 1.17 (55) and cf = 0.650 (49). Both 
data sets together can still be consistently fitted using 
the weighted average ci = 0.654 (49) of the leading order 
coefficients and 1-parameter fits for C2 in ([25| and ([26|) . 
This ensures that the ratio of the widths becomes one 
in the infinite volume limit. In contrast to the Polyakov 
loop maxima, the width of the transition is a physical 
observable, which is to leading order in the volume in- 
dependent of the BCs. We have chosen to show these 
I-parameter fits together with their Q values in Fig. [5] 



FIG. 6: Fits of the A*',- — 6 width of the transition. 



The methodology for the corresponding Nr ~ 6 fits 
is the same as before. Using the data of table IIVI we 



find the leading order coefficients 

'1 



1.27(11) and 

= 1.8(1.3), which average to ci = 1.27(11). With 
this values consistent 1-parameter fits are obtained and 
together with their Q values depicted in Fig. [6l Note 
that the ordinate scale in this figure is more than three 
times larger than in Fig. O Nevertheless the extracted 
physical values have to be the same to the extent that 
scaling holds. 

We want to plot the width in physical units of MeV 
versus the box size in fermi and follow a similar procedure 
as before for T'^(Ls). For given Ls{P^) we define 



AT^(L,)-T(/3^,- 



-A/32^/3/2)-r(/3p^,- 



-A/3,^/3/2), 



(27) 



The dependence AT{Ls) is shown in Fig. [T] Compared 
to periodic BCs disorder wall BCs lead to a substantial 
broadening of the transition for the volumes considered: 
At (lOfermi)^ by a factor of 4.3 and at (5fermi)'^ by 
a factor of 5.5. The width is slightly less than 1% of 
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N^=6 disorder wall BCs 
N^=4 disorder wall BCs 
N^=6 periodic BCs 
4 periodic BCs 
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11 



FIG. 7: Estimate of finite volume correction to the width of 
the deconfinement phase transition. 



the (enhanced) transition temperature at (lOfermi)^ and 
about 8% at (5 fermi)'^. Within the error bars, which are 
quite large for the widths, scahng works again weU. 

Simulations with disorder wall BCs have turned out to 
be far more CPU time consuming than those with peri- 
odic BCs. The decreased heights of our pseudo-transition 
signals (the maxima of the Polyakov loop susceptibili- 
ties), their increased widths and the strong finite size 
effects are the underlying reasons. While the reweighting 
range [H, [l^ is about the same for simulations with pe- 
riodic or disorder wall BCs on identically sized lattices, 
accurate disorder wall results require far more patches, 
i.e., independent simulations at distinct values. In ad- 
dition the signal is worse due to the decreased heights of 
the peaks. Finally extrapolations to infinite lattices, as 
needed for the finite volume continuum limit, are more 
demanding due to the strong and more sophisticated fi- 
nite size corrections. This is helped by including infinite 
volume extrapolations from periodic lattices as disorder 
wall data points, what can be done because these extrap- 
olations do not depend on the BCs. Still, data from our 
largest 6 x 48'^ lattice turn out to be essential to stabilize 
the iV^- = 6 fits for large lattice sizes. Far smaller lattices 
are sufficient when periodic BCs are used. 



tigation to a substantial rounding of the transition and to 
an increase of the (effective) deconfinement temperature 
by 5% to 20%. The physical reason for these corrections 
is that correlation lengths are proportional to and 
Lt I Ls does not approach zero in the finite volume con- 
tinuum limit. We estimate that the corrections are neg- 
ligible for the geometry of a torus, used in practically all 
previous LGT simulations of the subject, but not for BCs 
which reflect the cold exterior. Using the disorder wall 
BCs introduced in section III Al we find the magnitude 
of the effects on pure SU(3) LGT competitive to those 
of other corrections, foremost the inclusion of quarks. In 
particular the width of the transition increases by fac- 
tors 4 to 6 over the width found for a torus of the same 
physical size. 

Most scattering events at a RHIC are not from cen- 
tral collisions. So one has to cope with a distribution of 
volumes, each of it associated with its own effective de- 
confinement temperature and width. Due to such finite 
volume effects the concept of a sharp transition becomes 
blurred even when the effects of quarks, which convert 
the transition into a crossover [2(| , are not yet taken into 
account. Ultimately one may want to extend LGT stud- 
ies of the deconfining transition with BCs refiecting the 
confined outside world to full QCD. Before doing so, ad- 
ditional experience can be gained from pure SU(3) LGT. 
In future work we intend to perform simulation for the 
physically more realistic confinement wall BCs, which we 
introduced in section IIIBI Further, the entire equilib- 
rium thermodynamics of the dcconfined phase ought to 
be addressed. 

Close to the transition point the deconfined phase can 
only be studied non-perturbatively. Equilibrium QCD 
on the lattice allows this from first principles. But the 
assumption that equilibrium configurations can capture 
the essence of the transition is a strong one. Studies of 
the important dynamical aspects of the transition have 
to rely on phcnomcnological approaches. For instance, 
models which use hydrodynamics in the quark gluon 
plasma stage reproduce experimental data on particle 
abundances and flows well 1211. 
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